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Abstract 

We study mode-I fracture in lattices with noisy bonds. In contrast to previous attempts, by using 
a smah parameter that perturbs the force-law between the atoms in perfect lattices and using a 
3-body force law, simulations reproduce the qualitative behavior of the beyond steady-state cracks 
in the high velocity regime, including reasonable micro-branching. As far as the physical properties 
such as the structure factor g{r), the radial or angular distributions, these lattices share the physical 
properties of perfect lattices rather than that of an amorphous material (e.g., the continuous 
random network model). A clear transition can be seen between steady-state cracks, where a 
single crack propagates in the midline of the sample and the regime of unstable cracks, where 
micro-branches start to appear near the main crack, in line with previous experimental results. 
This is seen both in a honeycomb lattice and a fully hexagonal lattice. This model reproduces the 
main physical features of propagating cracks in brittle materials, including the behavior of velocity 
as a function of driving displacement and the increasing amplitude of oscillations of the electrical 
resistance. In addition, preliminary indications of power-law behavior of the micro-branch shapes 
can be seen, potentially reproducing one of the most intriguing experimental results of brittle 
fracture. 
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I. INTRODUCTION AND BRIEF REVIEW 



An extensive experimental effort on mode-I (tensile) fracture in amorphous materials 
has been made in the last two decades [llQ] (for a review, see [l3|). In the high crack 
velocity regime, the simple picture of a rapid steady-state crack generated via a given driving 
displacement, and exhibiting a given crack velocity (of order the Rayleigh surface wave 
speed) breaks down, and small micro-branches start to appear next to the main crack 
id . Il2| . Upon further Increase of the loading, the microscopic branches transform to large 
macro-branches . 

The experimental phenomena of the appearance of micro-branches has been the subject 
of an extensive theoretical effort as well. There have been several attempts within the 
framework of continuum models, based on the linear elasticity fracture mechanics (LEFM) 



theory 1^ . Yoffe predicted that steady-state cracks will become unstable at a specific crack 
velocity ~ 0.73c/j, based on maximal stress considerations 15|. Another ad hoc attempt 
based on energy considerations predicted a critical velocity of Vcr ~ O.Sci? [16]. However, 
the mode-I experiments have shown that the critical velocity is material-dependent, refuting 
those LEFM predictions. Several LEFM-based works predict a micro-branching instability at 
a material-dependent critical velocity \17], however, the specific parameter that determines 
the specific critical velocity, r(fcr) is an input parameter to the theory. 

Additional efforts have been made to explore the micro-branching phenomena based on 



LEFM 



18 



19| . However, as these works themselves argue, although they recover some 



features of the micro-branching instability, some main predictions such as that the energy 
flow to the micro-branches drops immediately, are unphysical (the micro-branch arrests im- 
mediately). These results raise the possibility that the micro-branches phenomena is a 3D- 
phenomena, while in 2D the high- velocity instability occurs only for extreme crack velocities 
{vcr > O.Scr), as the crack oscillates. This argument is based on experiments, both in biaxial 



mode-I crack in rubber 
ical works 
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20l |. and pure mode-I experiments in gels 2l|, and associated theoret- 



23| . However the mode-I experiments in PMMA seem to indicate that beyond 
a certain crack velocity, the crack structure that emerges is essentially two-dimensional 
There has been some success in exploring the instability using phenomenological mesoscale 



approaches based on the phase-field [2J] or on cohesive zones 
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27j : these however are 



difficult to quantitatively relate to an underlying microscopic picture. 



The failure of the continuum theory (LEFM) has given rise to an extensive theoretical 



effort using atomistic lattice models 
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44j . where 



and lattice simulations 

the inherent divergence of the elastic fields near the crack tip of the continuum theories is 
tamed, due to the finite atomistic lattice scale. These models yield steady state cracks 
without any additional parameters, once the inter-atomic forces are specified. Both lattice 
models and simulations have shown and reproduced the sharp transition between steady 
states cracks, where only the bonds on the midline of the lattice are broken, and the post- 
instability behavior, where beyond some critical velocity, other bonds start to fail. This 
critical velocity was found to have a strong dependence {0.3cr < Vcr ^ cr) on the parameters 
of the potential, i.e. it is material dependent. Although the lattice models yield the desired 
existence of a critical velocity, regarding the post-instability point behavior, the success is 



less impressive. Mode-III (out-o: 



patterns of micro-branches 
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?^lane shear mode) simulations have shown nice qualitative 
35l |. similar mode-I simulations (the mode for which most 



of the experiments actually have been performed), have failed to reproduce the qualitative 



39|. 



patterns of micro-branches emerging near the main crack [37, 

Thus, several attempts have been made to try to simulate cracks in amorphous materials, 
in which the bulk of experiments have been performed (e.g., glass, PMMA, Homalite-100 
etc.). The first attempts used the classic binary-alloy model (using two different kinds of 



particles) 45| simulating an amorphous material failed 46l-l48| : The crack always arrested. 
Recently 49| . a new approach were presented based on a continuous random network model 
(CRN) for simulating amorphous materials 50|, |5l|- In this model the sample looks like a 



distorted lattice, while each atom shares the same number of nearest neighbor atoms. The 
CRN model yields both steady-state cracks and the main features of the micro-branching 
instability, including the increasing size of the micro-branches and the increasing oscillations 
in the electrical resistance of the sample 491] with increasing external loading. Recently 45| . 
Dauchot, et al. succeeded in generating propagating steady-state cracks using the binary- 
alloy model, by going to the extreme brittle limit, where the force falls rapidly to zero for very 
small strains. However, no information was reported regarding the high-velocity instability. 
Moreover, the g{r) (RDF) generated from the amorphous model presented in 45|] are much 
less similar to the (?(r)'s of real amorphous materials than those the CRN's model generates. 

In this work we focus on trying to reproduce the successful results of the CRN model 
also in lattices (where the previous attempts have failed, as explained above). To do so, we 



let the force-law between the atoms vary slightly by changing randomly the lattice scale a 
between the atoms. As a result, the equilibrium locations of the atoms are slightly changed 
from their pure lattice locations, according to the modified force law. By breaking the pure 
symmetry of the perfect lattice, we hope to obtain a realistic micro-branching phenomena in 
these simple structures. The semi-quantitative behaviors of the micro-branches, such as the 
length of the micro-branches, should be less noisy than the corresponding results for CRN's. 

II. THE MODEL AND MAIN METHODOLOGY 

In our model, each bond (between atoms i and j) has a specific characteristic equilibrium 
distance at which the force is zero. This distance a^j is taken to vary slightly from the 
constant distance ao by a factor of e^j which is drawn from a uniform distribution: 

Uij ^ {1 + eij)ao, i^l,2,...,ne,toms,j&nn{i) (1) 

where tij G [—6, b] and 6 is a constant for a given lattice and in this work ranges between 
0^6^ 0.1, ttQ — 4 and nn{i) refers to the nearest-neighbors of site i. 

Between each two atoms there is a piece-wise hnear radial force (2-body force law) of the 
form: 

fl^j = kR{\fij\ - aij)9H (e - \fij\) fj^i, (2) 

where the Heaviside step function thetan guarantees that the force drops immediately to 
zero when the distance between two atoms \ fij\ reaches a certain value e > aij (the break 
of a bond). In this work we set £ = Oq -l- 1, and the units are chosen so that the spring 
constant kR is unity. Potentially, in addition there is a 3-body force law that depends on 
the cosine of each of the angles, defined of course by: 

— * — * 

\'i,j\\ii,k\ 

In a honeycomb lattice there are three angles associated with each atom and in the hexagonal 
lattice there are six of them (we note that in the hexagonal lattice this choice is a little bit 
arbitrary since there are in general additional optional angles for each atom). There is a 
certain preferred angle 9c for which the 3-body force law vanishes (in the honeycomb lattice 
we set 9c — and in the hexagonal lattice we set 9c — V^)- The 3-body force law that 
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acts on the central atom (atom i) of each angle may expressed as: 



fiXm = keicosOij^k - cosOc) 7^ Oh [e - \rij) n 



kg{cOs9ij k — COS 9c 



dfi 



(4) 



I /VI , . I I 'T' . 7 I 3 

\i i,k\ 



Oh {£ - \rij\) , 



while the force that is applied on the other two atoms (atoms j, k) may expressed as: 

Jj,{i,k) = ke{cos9ij^k - cosOc) - \ri,j\) r 



ke{cos6ij^k - cos 6c) 



(5) 



_ I ' i,J I I ' l,k I I I I I 



Oh {e-\ri,j\) 



Of course, the forces satisfy the relation: f^j^^^ = -{fj^{i^k) + fk,{i,j))- 

In addition, we used a Kelvin-type viscoelastic force, proportional to the relative velocity 



between the two atoms of the bond Vi^j: 



(6) 



with 7] the viscosity parameter. The viscous force vanishes after the bond is broken. Thus, 
the equation of motion of each atom is: 



j&pnn j,k(^3pnn ja&pnn 

where p = 1 for the honeycomb lattice and p = 2 for the hexagonal lattice (nn=nearest 
neighbors). The masses rrii can also be set to unity without loss of generality. 

The main methodology is as follows: After choosing the random value a^j for each bond, 
we allow the network to relax through a simple molecular-dynamics Euler scheme, in accord 
with Eqs. ([7]), with a non-zero rj, until the total energy is minimized. In Fig. [H we can 
see an example of a perturbed honeycomb lattice with b = 10%, while in Fig. [2] we can see 
the distribution of the radial distances of the bonds for a perturbed lattice with b = 2.5%, 
b = 10% along with the CRN, taken from 

We can see the qualitative difference between the perturbed lattice mesh and the CRN. 
While the radial distributions for the perturbed lattice mesh are flat in the range of 1 ± 6 
and then drops immediately to zero (since the random distribution was taken to be fiat in 



FIG. 1. A typical steady-state perturbed lattice grid using b = 10%. 



the range of 1 ± 6), the CRN has a long tails extending over larger distances. In addition, 
the angular distributions are much narrower than the CRN's angular distribution. 

A powerful tool to check the character of the grid is of course the radial distribution 
function (RDF or g{r)). In Fig. [3] we can see the RDFs of the perturbed lattices and 
the CRN. We can see again the qualitative difference between the meshes. While the CRN 
looks very much like a real amorphous material (see 



the RDFs of the perturbed lattices 
look exactly like a pure lattice RDF (set of (5-functions) , only slightly perturbed due to the 
random noise, even for large r. Thus, the structure of the perturbed lattices is like a lattice 
material rather than an amorphous material. 

After characterizing the initial lattices, we stretch the lattice under a mode-I tensile 
loading with a given constant strain using a given driving displacement A, and seed the 
system with an initial crack. We let the crack propagate via the same molecular dynamics 
Euler scheme that was introduced before. The lattice mesh we use contains 162-272 ^ 45, 000 
(A^ = 80 in the Slepyan model notation) atoms for the honeycomb lattice and 162 ■ 408 ~ 
65, 000 atoms for the hexagonal lattice. 



III. HONEYCOMB LATTICE 

Using 6 = 0, i.e. an unperturbed honeycomb lattice, we obtain the well-known non- 
physical behavior of the crack above threshold. For small strains we get a perfect steady- 
state crack, while upon increasing the driving displacement, the crack bifurcates to two 



macro-branches that propagate to the edges of the sample 



371]. The same happens in the 
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FIG. 2. The radial (a) and the angular (b) distributions for a perturbed lattice mesh with b = 2.5%. 
(c-d) Same for b = 10%. (e-f) Same for the CRN amorphous mesh. 



honeycomb lattice including the 3-body force law, both with large-viscosity and with a 
negligible viscosity (Fig. H]). 

Using a finite value of b to perturb the lattice, we obtain nice snapshots of micro-branches, 
very much like those obtained using the CRN 49|. In Fig. |5]^a) we can see that when the 
driving displacement exceeds some value, a large micro-branch starts to appear. In Fig. [6] 




FIG. 3. The radial distribution function g{r) as a function of r (in lattice scale units) in the 
perturbed-lattice systems with different magnitude and in the continuous random network system. 
A major qualitative difference can be seen between the two networks. 




FIG. 4. The propagating crack snapshot using a pure honeycomb lattice using ij = 0.25 and 
'^/Ag = 2.7. The crack bifurcates to two that travel to the end of the sample. 

we can see the final pattern of broken bonds for two cases, one for t] = 2 (Fig. [6]^a)) and 
one for rj = 2.5 (Fig. [6]^b)). The patterns looks very much like the fracture pattern seen 
using the CRN. Moreover, the micro-branches also look similar to the experimental images 



of micro-branches in PMMA js-llOl]. This is an important result. As far as we know, this is 
the first time that such a micro-branch pattern has appeared in a lattice material (previously 
shown by us only using an amorphous material model 49(]). We note that when the main 
crack continues and the micro-branch arrests, one piece of the lattice overlaps with another 
piece (see Fig. Et^b)). This is a non-physical effect and is caused from the fact that cracking 
in this model is irreversible. This effect is not large using the honeycomb lattice, but will 
be much more pronounced for the hexagonal lattice. 

Beyond obtaining the qualitative features of the micro-branches, the perturbed lattice 
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FIG. 5. (a) A propagating crack snapshot using a perturbed honeycomb lattice model using 
b = 10%, 7] = 2 and ^/Ag = 3.8. (b) A snapshot of the overlapping problem. After cracking, the 
two pieces of the cracked lattice overlap, because fracture is an irreversible process in this model. 




1 




FIG. 6. (a) The pattern of the final cracked lattice using a perturbed honeycomb lattice model 
using b = 10%, r/ = 2 and ^/Aq = 3.8. (b) Same with r/ = 2.5 and ^/Aq = 4. 



model reproduces the main semi-quantitative results of mode-I fracture beyond the onset of 
instability. Because of the chaotic nature of the problem (small changes in the simulation 
parameters yield different crack patterns, but with similar quantitative properties, such 
as the crack velocity), we changed the time-step by a little bit (±15%) for each driving 
displacement A, to have sufficiently good statistics on the resulting parameters (about 20 
runs for each point, therefore each point in Fig. [71^a-c) represents a set of ~ 20 runs). In 
Fig. [7] we can see the v{A) curve (a), the total length of the micro-branches (b) and the 
amplitude of the oscillations of the electrical resistance (c). The v(A) curve looks very much 
like the typical v{A) curve for mode-I fracture (for example, see The error 

bars represent the statistical error of the crack's velocity using several simulations for each 
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FIG. 7. (a) The resulting v{A) curve for a perturbed honeycomb lattice model using b = 10% with 
different rj. The error bars were calculated using several simulations, each with small change in dt, 
due to the chaotic nature of the problem (b) Size of total number of micro-branches as a function 
of the crack velocity v. For small velocities the total size of micro-branches tends to zero (a steady 
state crack), (c) The amplitude of the oscillations in the electrical resistance as a function of the 
crack velocity. Even in steady-state cracks, the amplitude of the oscillations tends to a certain 
finite (non-zero) value. 



A. The graph of the total length of the micro-branches as a function of the crack's velocity 
is one of the most important results of this work. We can easily see that the total length 
of micro-branches goes to zero at small velocities, with a clear growth (near v ~ 0.57cr) 
with increasing velocity (or A). This is in direct accord with the to experiment results (but 
instead of a sharp transition between steady-state cracks area and micro-branches area, 
we get a smooth transition, due to the noisy character of discrete atomistic simulation at 
smaller scales). Actually, the experimental results refer as the length of an average single 
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micro-branch, but because the lack of statistics (we have only a few micro-branches in each 
single simulation, so it is hard to define the length of an average single micro-branch), we use 
the total sum of broken bonds instead (besides of course the main crack). This should be a 
sufficiently close substitute (we see that we do noi have just more "short" micro-branches at 
large drivings; the size of each micro-branch indeed grows). Comparing to the corresponding 
result using the CRN for amorphous material 49(| (see also the triangles in Fig. |8]^b)), where 
this transition was less pronounced, here it is much more clear. We note that by defining 
the total size of the micro-branches we subtracted all the "micro-branches" of size 1 or 2 
broken bonds, which we neglect and treat as numerical noise. 

From Fig. El^b) we can see that using t] = l.b yields a too noisy system and we get either 
arrested cracks or significant micro-branches. There is no intermediate zone of steady-state 
cracks. Only using r/ = 2 do we yield clear steady-state cracks. Increasing rj further more to 
?7 = 2.5 does not change the results appreciably. In addition, the amplitude of the oscillations 
of the electrical resistance shows a nice agreement with the experimental result, as well as 
the CRN results 49| (For a wide discussion in the different terminology between the crack's 



velocity oscillation and the electrical resistance oscillations and the appropriateness of the 



"electrical resistance" oscillations as a diagnostic, see j49|). Beyond the "critical velocity" 
the amplitude of the oscillations increases rapidly, while for small velocity the amplitude of 
the oscillations is constant. 

The sensitivity to the value of 6, characterizing the width of the bond length distribution, 
was also investigated, exploring how much we can reduce h and still get a physical behavior 
of micro-branches pattern, recalling that 6 = (pure lattice) does not create a physically 
realistic pattern. For h ^ 1%, we reproduce the non-physical perfect lattice behavior, i.e. 
the perfect lattice behavior is not a singular case, in the sense that infinitesimal change in 
the lattice produce micro-branches. You need a significant perturbation to yield physical 
behavior. In fact, even for h = 2.5%, the micro-branches seems to be in an almost straight 
lines, along the preferred lattice directions, and thus, less physical. In Fig. [8] we can see the 



quantitative results using different values of 6, along with the CRN results (taken from |49| ) . 



From the f (A) curve (Fig. El^a)) and the amplitude of the oscillations (Fig. El^c)) we can 
see that the CRN is closer to the b = 10% results. In Fig. [8]^b) we see that low b results look 
like the CRN result. This fact encourages us to conclude that the transition of the CRN 
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FIG. 8. (a) The resulting w(A) curve for a perturbed honeycomb lattice model with different values 



491]) using rj = 2. The CRN curve is closer to the 



of b, along with the CRN results (taken from 
b = 10% curve (b) Size of total number of micro-branches as a function of the crack velocity v. 
The somewhat questionable transition between steady-state cracks and micro-branches behavior 
using the CRN, looks much more clear using perturbed lattices models, (c) The amplitude of the 
oscillations in the electrical resistance as a function of the crack velocity. Again, the CRN curve is 
closer to the b = 10% curve. 



results between low velocities and high is real, since increasing b we get the same effect, but 
with a much sharper transition. 



IV. HEXAGONAL LATTICE 



The classic models concernin g rn olecula r dy namics fracture simulations in perfect lattices 
used a hexagonal lattice (30|, ImI, l36M39l l4ll |43|, |4J]. As mentioned in the Introduction, those 
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FIG. 9. The radial (a) and the angular (b) distributions for a perturbed hexagonal lattice mesh 
with and without 3-body force law using different values of b. kg/kr = 1 and b = 10% yields the 
almost the same distributions as kg/kr = and b = 10% does. 



models were able to show nice qualitative results in the mode-Ill fracture simulations 30l.l35|. 
but failed to create the physical pattern of micro-branches in mode-I fracture simulations, 
investigated in the experiments 37|, |39|. The new findings from the previous section, that a 



small perturbation in the potential between each two atoms can create a physical pattern of 
mode-I fracture, suggest it is worthwhile to try it also in the classic hexagonal lattice. The 
resulting radial and angular distributions for several parameters are shown in Fig. |9l 

Using the same magnitude of perturbation as in the honeycomb lattice, b = 10% (with 
no 3-body force law, i.e. kg = 0) yields an extreme noisy simulations, and eventually, in 
most cases, the crack arrests. Using a smaller value of 6 (6 = 2.5%), yields in general micro- 
branching patterns (for large driving displacement) (Figs. ITO] and ITTl) . but it still looks less 
physical than the honeycomb lattice results. At low drivings, there is a single steady-state 
crack propagating in the midline of the sample yielding no broken bonds besides the main 
crack. Nevertheless, the benefit of using this model (in the absence of a 3-body force law) 
is that we can compare it to Slepyan's lattice steady-state models (like in 3^, using large 
a, which corresponds to a piecewise-linear model), when the origin of instability between 
steady-state cracks and micro-branches behavior is known exactly (although the models use 
6 = 0, the small perturbation does not change the results significantly). Using b smaller than 
2.5% reproduces the non-physical behavior of perfect lattices. In addition, the hexagonal 



mesh allows us to work with small values of rj which is more relevant experimentally 53l. |54|. 
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FIG. 10. (a) A propagating crack snapshot using a perturbed hexagonal-lattice model using b = 
2.5%, Tj = 0.25 and ^/Aq = 2.2. (b) A snapshot of an extensive overlapping problem. After 
cracking, the two pieces of the cracked lattice overlap, because fracture is irreversible process in 
this model. 




FIG. 11. The pattern of the final cracked lattice using a perturbed hexagonal-lattice model using 
b = 2.5%, r? = 0.25 and a/aq = 2.2. 



than the large 1] honeycomb lattice. 

In Fig. [TOT a) we can see the a birth of a micro-branch that bifurcates from the main crack, 
that eventually arrests. The pattern looks very much alike the honeycomb lattice. In Fig. 
[11] we can see the final pattern of broken bonds for the case of rj = 0.25. The quantitative 
results regarding the f (A) curve and the total size of micro-branches as a function of v are 
presented in Fig. [121 

We can see in Fig . IT2\ &) that the velocities below v = cr reproduce the Slepyan's lattice 
model results from [37|, yielding perfect steady-state cracks, with no micro-branches at 
all (Fig. [T2r b)). Increasing the driving displacement A further, yields a non steady-state 
behavior, as the steady-state lattice model solution becomes unstable, yielding a micro- 
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FIG. 12. (a) The resulting w(A) curve for a perturbed hexagonal-lattice model using different 
values of kg with r] = 0.25, (6 = 2.5% for ke/kr = 1 and 6 = 10% for feg/A;^ = 10.5 = 10%. 
The error bars were calculated using several simulations, each with small change in dt, due to the 
random nature of the problem (b) Size of total number of micro-branches as a function of the 
crack's velocity v. In small velocities the total size of micro branches tend to zero (a steady state 
crack). The including 3-body results are normalized to the non 3-body results. 



branching behavior (again, only above a threshold h). We note again that in many cases 
the fracture pattern do not looks as physical as in Fig. [11], yielding some non-physical 
results. In addition, the problem of overlapping zones is more extensive in the hexagonal 
lattice (see Fig. [TOlfb)). yielding large areas of overlapping zones and thus, yielding an 
unphysical behavior. The main benefit so far of using this model is that Fig. [T2](b) yields 
similar qualitative results to Fig. El^b) or Fig. [8]^b), emphasizing that the transition at low 
velocities to a steady-state behavior is real, since in a hexagonal case we get absolutely zero 
micro-branches at low-velocities, and agreement with the Slepyan's lattice models. 

Adding a 3-body force law in a hexagonal lattice, one must use large values of kg since 
the relatively large number of nearest neighbors doesn't allow each 3 atoms to generate an 
angle significantly larger than ^/-i using h/Q ~ (rs was the case in the honeycomb lattice), 
and thus, the 3-body energy is negligible. When we increase kg using h = 2.5%, we yield 
a perfect lattice behavior because kg is too large. Increasing b further to 6 = 10% yields a 
good balance which the 3-body energy is not negligible and yields similar radial and angular 
distributions as a hexagonal lattice without the 3-body force law using b = 2.5% (see Fig. 
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FIG. 13. The crack pattern for a perturbed hexagonal-lattice model using kg/k^ = 10, b = 10% 
and rj = 0.25. (a) ^/Aq = 1.8 with different vahies of dt. (b) ^/Aq = 1.9 with different values of dt 
for the first two upper patterns, and ^/Ag = 2 for the bottom pattern. 



We normalize the f (A) curve results and the total size of micro-branches as a function of v 
using this model to the values of the non-3-body force law case (the normalization constants 
Ac and cr are of course larger using large significant kg) and present them in Fig. [T21 We 
can see that these two models share the same qualitative behavior (although the including 
3-body force law results yields a small number of micro-branches at low velocities, very 
much like the honeycomb case). In addition, this model still suffers from a severe problem of 
overlapping of pieces of the mesh after branching (in contrast to honeycomb lattice, which 
in this model, this problem is minor). 

However, two surprising results appear using this model (a perturbed hexagonal lattice 
including 3-body force law). First, the main crack stays more confined to the middle of the 
sample, even for large driving, with large micro-branches, in contrast to all other models, 
including the honeycomb lattice model. Second, and most important, the larger micro- 
branches have a non-linear, power law shape (!) very much alike the experimental data. In 
Fig. [T3]we can see several of crack-patterns for different driving displacements (with different 
dt, with shifting in the y-axis, of course). We add a basic power law fit (with different xq 
and Uq for each large micro-branch, of course). We got a nice power law behavior (but with 
power of 0.5 — 0.65 instead of ~ 0.7 of js, 7|). This is of course a very preliminary result, 
and must be tested in larger scales. We note that the ~ 0.7 power law behavior is not 



universal for all the experiments. At least in the experiments of [l^ (Fig- 3), 55| (Fig. 
4(a)) and 12|] (Fig. 1(b)), straight micro-branches appear, so the physical behavior of Figs. 
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[6] and [TT] is important too. On the other hand, power law behavior (different than 1) is 
seen in the atomistic model only in a perturbed hexagonal lattice including 3-body force 
law. These observations support the conclusion that the macroscopic behavior of fracture 
depends strongly on the inter-atomic microscopic potential. 

V. DISCUSSION 

We have shown that micro-branching can be reproduced in lattice materials, using a 
small perturbation parameter b which perturbs the inter-atomic potential between each two 
atoms. In addition to the qualitative patterns of micro-branches, semi quantitative results 
are shown, particular the total size of micro-branches (which corresponds to the average 
size of a micro-branch in the experiments) as a function of the velocity. A clear transition 
between steady-state behavior and the post-instability region is seen, characterized by an 
increased number of broken bond off the midline of the sample. This result is in line with 
the results of the CRN model. In particular we obtain in the hexagonal case (without 3- 
body force law) no micro-branching at all, as in the Slepyan lattice models. The increased 
amplitude of the RMS of the electrical resistance is shown as well, in agreement to the 
experiments. 

In addition, preliminary signs of power law behavior for the shape of the side branches can 
be seen using a hexagonal lattice including a 3-body force law, in accord with experiment, 
and which has not been observed in previous atomistic models. Much more extensive work 
is underway to extend this result to larger scales, as well as exploring the opening angle of 
the crack. 
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